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There  is  a  significant  push  to  increase  the  output  power  of  thermoelectric  generators  (TEGs)  in  order  to 
make  them  more  competitive  energy  harvesters.  The  thermal  coupling  of  TEGs  has  a  major  impact  on  the 
effective  temperature  gradient  across  the  generator  and  therefore  the  power  output  achieved.  The 
application  of  micro  fluidic  heat  transfer  systems  (pHTS)  can  significantly  reduce  the  thermal  contact 
resistance  and  thus  enhance  the  TEG’s  performance.  This  paper  reports  on  the  characterization  and 
optimization  of  a  pTEG  integrated  with  a  two  layer  pHTS.  The  main  advantage  of  the  presented  system  is 
the  combination  of  very  low  heat  transfer  resistances  with  small  pumping  powers  in  a  compact  volume. 
The  influence  of  the  most  relevant  system  parameters,  i.e.  microchannel  width,  applied  flow  rate  and  the 
pTEG  thickness  on  the  system’s  net  output  performance  are  investigated.  The  dimensions  of  the  pHTS/ 
pTEG  system  can  be  optimized  for  specific  temperature  application  ranges,  and  the  maximum  net  power 
can  be  tracked  by  adjusting  the  heat  transfer  resistance  during  operation.  A  system  net  output  power  of 
126  mW/cm2  was  achieved  with  a  module  ZT  of  0.1  at  a  fluid  flow  rate  of  0.07  1/min  and  an  applied 
temperature  difference  of  95IC 

It  was  concluded  that  for  systems  with  good  thermal  coupling,  the  thermoelectric  material  optimi¬ 
zation  should  focus  more  on  the  power  factor  than  on  the  figure  of  merit  ZT  itself,  since  the  influence  of 
the  thermal  resistance  of  the  TE  material  is  negligible. 

©  2013  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Due  to  the  simultaneous  rise  in  energy  consumption  and  envi¬ 
ronmental  awareness,  the  worldwide  demand  for  more  efficient  and 
clean  energy  systems  is  growing.  One  promising  approach  to 
improve  a  system’s  efficiency  is  to  recover  the  produced  waste  heat 
by  means  of  the  thermoelectric  effect.  Thermoelectric  power  gen¬ 
eration  can  potentially  be  applied  for  waste  heat  recovery  in  energy 
conversion  systems  [1],  industrial  processes  [2]  or  automotive  ap¬ 
plications  [2,3].  So  far,  however,  the  commercial  applications  of 
thermoelectric  generators  (TEGs)  have  been  limited  to  niche  mar¬ 
kets,  such  as  space  [4,5]  and  remote  or  hazardous  places  [6,7],  The 
main  reasons  for  this  are  (1)  the  relatively  low  thermoelectric  con¬ 
version  efficiency,  (2)  high  fabrication  costs  of  thermoelectric  mod¬ 
ules,  and  (3)  suboptimal  exploitation  of  the  available  temperature 
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gradients.  While  much  effort  is  spent  to  improve  the  thermoelectric 
figure  of  merit  (ZT)  by  quantum  confinement  or  phonon  scattering 
[8—10]  as  well  as  to  develop  low-cost  fabrication  technologies  [11], 
less  attention  has  been  paid  to  the  optimization  of  the  overall  sys¬ 
tem’s  performance  so  far.  This  includes  thermal  coupling  of  the 
thermoelectric  device  to  the  heat  source  and  heat  sink,  as  well  as  an 
optimal  matching  of  the  thermal  contact  and  TEG  resistances  [  12],  By 
decreasing  the  thermal  contact  resistance  to  the  cold  and  hot  side 
reservoir,  a  significant  enhancement  of  the  generator’s  output  power 
can  be  achieved  [13],  The  major  impact  of  the  thermal  contact 
resistance  on  the  output  performance  of  a  thermoelectric  generator 
is  illustrated  in  Fig.  1  and  compared  to  the  effect  of  ZT.  While  the 
influence  of  ZT  on  the  output  power  has  a  close  to  linear  dependency, 
the  thermal  contact  resistance  exhibits  a  power  dependency  pro¬ 
portional  to  1  /R2on. 

The  most  straight-forward  approach  to  achieve  small  thermal 
contact  resistances  is  to  apply  fluidic  heat  exchanger  systems  for 
the  heat  supply  and  dissipation.  This  is  also  the  most  effective 
method  for  waste  heat  recovery.  Several  large  scale  fluidic  systems 
using  liquid  as  active  media  have  been  tested  for  use  in  industrial 
applications.  Tsuyoshi  et  al.  [14]  reported  experimental  results  on  a 
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Fig.  1.  Simulated  output  power  of  a  215  pm  thick  Bi2Te3  pTEG  (RTf.g  =  2.8  cm2K/W)  in 
dependence  of  the  ZT-value  and  the  thermal  contact  resistance  Rcon  at  an  applied 
temperature  difference  of  20K.  The  simulations  were  performed  applying  a  thermo¬ 
electric  model  under  matched  electrical  load  [17], 

thermoelectric  engine  made  out  of  TEGs  stacked  between  a  parallel 
plate  heat  exchanger  using  oil  and  water  as  active  media.  At  an 
applied  temperature  gradient  of  130K,  an  output  power  of  170  W 
could  be  reached.  Niu  et  al.  [15]  built  a  similar  parallel  plate  heat 
exchanger  with  commercially  available  Bi2Te3  modules  reaching, 
140  W  at  an  inlet  temperature  difference  of  120K.  Crane  et  al.  [2] 
constructed  a  TEG-heat  exchange  assembly  with  a  stack  of  6  TEG 
modules  producing  500  W  at  an  inlet  temperature  difference  of 
205K.  Due  to  the  growing  interest  in  thermoelectric  waste  heat 
recovery  for  automotive  applications,  several  gas/liquid  heat  ex¬ 
change  systems  have  been  studied  and  tested  [3].  First  functional 
prototypes  have  already  been  installed  in  test  engines  and  values  as 
high  as  507  W  at  a  AT  of  553K  could  be  reached  [16], 

Although  the  achieved  output  power  values  from  the  experi¬ 
mental  studies  are  remarkable,  the  systems  are  generally  very  large 
and  heavy  resulting  in  rather  low  output  powers  per  volume.  Often, 
standard  large  scale  heat  exchangers  are  applied,  where  high 
amounts  of  the  working  fluid  are  pumped  through  the  system  and 
the  consumed  pumping  power  is  neglected.  Furthermore,  the  re¬ 
ported  systems  usually  work  without  optimized  thermal  contact 
resistance  and  overall  system  performance.  In  contrast,  scaling  down 
the  heat  exchangers  to  the  micrometer  range  brings  several  advan¬ 
tages.  Micro  heat  transfer  systems  (pHTS)  can  achieve  high  heat 
fluxes  due  to  an  increase  of  the  convective  heat  transfer  coefficient  as 
well  as  of  the  surface  area  [17],  Even  though  pHTS  generally  exhibit 
larger  hydrodynamic  resistances,  the  effect  on  the  pumping  power 
can  be  compensated  by  smart  design  and  smaller  flow  rates  needed. 
With  the  reduced  heat  transfer  resistances,  thinner  generators  can  be 
applied  for  thermal  resistance  matching.  This  makes  the  system 
compatible  for  waste  heat  recovery  in  industrial  applications,  where 
high  heat  fluxes  are  often  necessary.  Additionally,  more  compact 
systems  with  reduced  size  and  weight  can  be  built,  resulting  in 
increased  applicability  and  modularity  of  the  system.  The  perfor¬ 
mance  characteristic  of  different  micro  heat  sinks  for  chip  or  LED 
cooling  applications  have  been  investigated  in  several  theoretical  and 
experimental  studies  [17,18],  Rezania  et  al.  [19]  evaluated  the  theo¬ 
retical  performance  of  a  thermoelectric  generator  using  a  rectangular 
microchannel  heat  sink.  The  simulation  results  showed  no  net  TEG 
power  improvements  of  the  micro  heat  sink  system  compared  to 
macro  scale  heat  sinks,  mainly  due  to  the  large  pressure  drops  inside 
the  microchannels.  To  avoid  such  large  pressure  losses  and  to  still 
benefit  from  the  improved  thermal  performance  of  microchannels, 
multiple  layer  manifold  systems  can  be  applied  [20—22], 


Fig.  2.  Exploded  view  of  the  pHTS/pTEG  system. 

Therefore,  a  novel  approach  of  minimizing  thermal  resistances  by 
integrating  two  layer  pHTS  with  pTEGs  is  proposed.  This  paper 
presents  an  experimental  study  on  the  performance  of  such  an  in¬ 
tegrated  system  as  shown  in  Fig.  2.  The  pHTS  consists  of  a  copper 
microchannel  layer  for  efficient  heat  removal  and  a  polymeric 
manifold  channel  layer  responsible  for  the  fluid  distribution  [13]. 
Besides  a  small  pressure  loss  and  a  low  heat  transfer  resistance,  the 
proposed  system  also  benefits  from  uniform  cooling  properties. 
Many  different  sets  of  parameters  influence  the  performance  of  the 
described  pHTS/pTEG  system;  namely  the  dimensions  of  the  pHTS 
itself,  the  properties  and  dimensions  of  the  thermoelectric  generator 
and  the  applied  boundary  or  working  conditions.  The  paper  focuses 
on  the  analysis  and  characterization  of  the  most  important  pHTS/ 
pTEG  system  parameters  with  the  aim  of  optimizing  the  overall 
system’s  net  output  performance.  In  the  first  part,  the  pHTS  is 
characterized  in  regards  to  the  system’s  dimensions,  flow  rates  and 
applied  temperature  gradients.  In  the  second  part,  the  output  power 
of  the  pHTS/pTEG  systems  is  investigated  with  respect  to  the  pHTS 
dimensions,  boundary  conditions  and  pTEG  thermal  properties.  The 
advantages  of  variable  thermal  resistances  are  discussed,  and  prac¬ 
tical  aspects  of  the  thermal  resistance  matching  are  evaluated. 

2.  Experimental 

The  pHTSs  were  designed  based  on  system  and  fabrication 
constraints  in  combination  with  performed  simulations  results.  The 
parameter  with  the  largest  impact  on  the  thermal  and  hydrody¬ 
namic  performance  of  the  pHTS  was  identified  to  be  the  micro- 
channel  dimensions.  Therefore,  pHTSs  with  different  channel 
widths  were  fabricated  and  analyzed.  The  relevant  parameters  of 
the  pHTSs  used  are  summarized  the  Table  1. 

The  investigated  pHTSs  were  fabricated  using  standard  photo¬ 
lithography  and  electrochemical  deposition.  The  detailed  working 
principle  and  the  fabrication  process  of  the  pHTS  are  reported  in 
Ref.  [13],  The  pTEGs  used  were  provided  by  greenTEG  GmbH  and 
are  based  on  the  electrodeposition  of  Bi2Te3  thermoelectric  ma¬ 
terial  into  flexible  polymer  molds.  A  detailed  fabrication  process  is 

Table  1 

pHTS  design  parameters. 

MicroChannel  MicroChannel  Fin/channel  Manifold  Manifold  Total 

width  -  wmch  height  -  hmch  ratio  [-]  channels  channels  area 

[pm]  [pm]  width  -  wmf  height  -  hmf  [mm] 

h™>l [pm] 

35, 50, 60, 80  190  1  150-300  1000  8  x8 

(tapered) 
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described  in  Ref.  [11].  In  order  to  provide  electrical  isolation,  a 
2.5  pm  layer  of  photoresist  (Microposit  S1813,  Shipley)  was 
deposited  as  a  passivation  layer  on  each  side  of  the  pTEG.  To 
minimize  the  interface  resistance  to  the  pHTS  and  the  heater,  the 
passivated  generator  was  additionally  covered  with  approximately 
30—40  pm  thermal  paste  (Dow  Corning  TC-5026).  Prior  to  their 
integration  into  the  pHTS,  the  pTEG-units  were  characterized  in  a 
thermoelectric  measurement  setup  [23].  The  relevant  parameters 
of  the  standard  pTEG  used  are  summarized  the  Table  2.  The 
module  ZT  is  defined  as  the  figure  of  merit  of  the  generator,  where 
the  complete  pTEG  resistance  (including  electric  contact  resistance 
and  the  interconnect  resistance)  is  used  for  the  electric  resistivity 
calculation. 

For  the  characterization  of  pHTSs  and  pHTS/pTEG  systems,  the 
devices  were  mounted  on  a  clamping  fixture  and  connected  to  a 
controlled  fluid  loop  [13].  The  clamped  pHTS/pTEG  system  and  its 
equivalent  thermal  circuit  model  are  schematically  depicted  in 
Fig.  3.  The  measurement  fixture  contains  a  resistive  heater  (plat¬ 
inum  meander  on  AI2O3  substrate),  a  copper  block  with  an  inte¬ 
grated  temperature  sensor  and  a  pressure  controlled  clamping  tool. 
The  equivalent  thermal  circuit  of  the  pHTS/pTEG  systems  consists 
of  the  heat  transfer  resistance  of  the  pHTS  RHts.  the  thermal 
resistance  of  the  generator  Rteg  and  the  interface  resistances  Rin- 
terface  accounting  for  the  electrical  passivation  layer  of  the  TEG  and 
the  thermal  paste.  Therefore,  the  total  contact  resistance  is  defined 
as  the  sum  of  the  variable  Rhts.  and  the  two  ^interface- 

The  cooling  liquid  (DI  water)  flow  was  measured  by  means  of  a 
Coriolis  flow  meter  (error:  ±1%  of  reading).  The  pressure  drop 
across  the  pHTS  was  recorded  with  a  differential  pressure  sensor 
(error:  ±1  mbar).  The  consumed  pumping  power  was  calculated 
from  the  product  of  the  fluid  flow  and  the  measured  pressure  drop. 
The  heat  flux  through  the  system  was  defined  by  the  power 
consumed  by  the  temperature  controlled  heater.  The  heat  transfer 
resistance  of  the  pHTS  was  calculated  by  dividing  the  applied 
temperature  difference  AT  by  the  measured  heat  flux.  The  ther¬ 
moelectric  power  output  was  calculated  from  the  measured  See- 
beck  voltage  VSb  under  electrically  matched  loads  conditions  and 
corresponds  to  the  maximal  achievable  power  output  [24],  A  fluid 
inlet  temperature  of  293K,  a  tool  pressure  of  660  kPa,  a  measure¬ 
ment  duration  of  120  s  and  a  sampling  frequency  of  0.76  s-1  were 
defined  as  standard  test  conditions. 

Before  the  performance  measurements  of  the  pHTS/pTEGs  sys¬ 
tems,  the  pHTSs  were  characterized  with  respect  to  their  thermal 
and  hydrodynamic  performance  inside  the  fluid  loop  measurement 
setup.  To  investigate  optimal  system  dimensions,  pHTSs  with  four 
different  microchannel  widths  were  analyzed  (Table  1 ).  Addition¬ 
ally,  the  influence  of  external  boundary  conditions  on  the  pHTS 
performance,  i.e.  the  fluid  inlet  temperature  and  heater  tempera¬ 
ture,  was  studied.  The  net  output  performance  of  the  pHTS/pTEG 
system  was  characterized  with  respect  to  the  fluid  flow  rate,  the 
applied  temperature  gradient  and  the  microchannel  dimensions. 
The  influence  of  thermal  resistance  matching  by  pTEG  thickness 
variations  was  also  investigated.  For  this  purpose,  the  pTEG  was 
ground  down  to  three  different  heights  and  re-contacted  after  each 
grinding  step. 

The  measurement  uncertainty  was  calculated  according  to  ISO/ 
IEC 17025  standards  [25],  The  resulting  standard  deviation  from  the 


Measured  B&Tea  |iTEG  parameters. 


Module  ZT  Electrical 

Thermal 

Total 

Ratio 

[-]  resistance  Rei 

resistance  Rteg 

thickness 

passive/TE 

m 

[cm2K/W] 

[pm] 

material  [— ] 

0.06  ±  0.01  2.35  ±  IE-3 

5.1  ±  0.3 

215  ±10 

1.43  ±  0.1 

ll  l  Hill  Tqoj, 
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Fig.  3.  Schematic  of  the  pHTS/pTEG  system  clamped  inside  the  measurement  fixture 
and  its  equivalent  thermal  circuit  model. 


time  averaged  measurement  signals  (random  error)  was  combined 
with  measurement  equipment  errors  (systematic  errors).  The 
Gauss  propagation  of  uncertainty  was  used  for  error  calculations  of 
indirectly  measured  values  (Rhts.  Ppump.  Pteg).  The  measurement 
repeatability  of  the  clamping  structure  was  tested  by  multiple 
mountings  of  the  same  pHTSs  and  pHTS/pTEG  devices.  The 
maximum  measured  standard  deviation  errors  were  below  5.3% 
and  0.52%  for  the  measured  heat  transfer  resistance  Rhts  and  See- 
beck  voltage  Vsb.  respectively.  This  measured  standard  deviations  of 
the  clamping  influence  are  all  below  the  calculated  measurement 
uncertainty  plotted  as  error  bars  in  the  figures. 

The  pHTS  model  used  was  a  semi-empirical  one-dimensional 
model  [26]  extended  with  a  discrete  temperature  distribution 
along  the  fluid  path  inside  the  microchannels  and  temperature 
dependent  fluid  parameters.  The  modeled  Rhts  is  composed  of 
three  thermal  resistances  in  series:  the  conductive  resistance 
through  the  microchannels,  the  convective  resistance  at  the  fluid- 
channel  interface  and  the  fluid  resistance  due  to  the  limited  heat 
capacity.  The  resistance  of  the  thermal  paste  was  fitted  to 
0.08  cm2K/W  corresponding  to  a  layer  thickness  of  30  pm.  For  flows 
below  0.4  1/min,  a  correction  factor  was  introduced  in  order  to 
account  for  a  reduced  heat  transfer  due  to  randomly  trapped  air 
bubbles  in  the  microchannels  and  possible  heating  up  of  the  fluid 
inside  the  manifold  channels.  For  the  pHTS/pTEG  system,  the  ID 
model  was  extended  with  a  TEG  model  under  matched  electrical 
load,  taking  Joule  heating  and  Peltier  effects  into  account  [27], 
Interface  resistances  accounting  for  the  electrical  passivation  and 
the  thermally  conductive  paste  were  added  (0.22  cm2I</W  at  each 
side).  The  measured  pTEG  and  pHTS  parameters  from  Tables  1  and  2 
were  used  for  the  simulations. 


3.  Results  and  discussion 

3.1.  fiHTS  characterization 

The  thermal  and  hydrodynamic  performance  of  the  pHTS  with 
different  microchannel  dimensions  was  investigated,  and  the  in¬ 
fluences  of  varying  boundary  conditions  were  analyzed.  Fig.  4a) 
shows  the  measured  and  computed  heat  transfer  resistance  Rhts  of 
the  pHTSs  with  different  microchannel  dimensions  relative  to  the 
fluid  flow  rate.  The  heat  transfer  resistance  decreases  with  increasing 
flow  rates  mostly  due  to  a  rise  in  the  heat  capacity  rate  of  the  fluid. 
The  curves  flatten  down  towards  higher  flow  rates,  since  the  relative 
influence  of  an  increased  effective  mass  flow  on  the  heat  flux  is 
reduced.  With  decreasing  channel  widths,  Rhts  is  reduced  mainly  due 
to  an  increase  in  the  convective  heat  transfer.  At  smaller  micro- 
channel  widths,  the  thermal  boundary  layer  is  thinner  and  therefore 
the  heat  convection  inside  the  channel  is  improved.  Additionally,  the 
effective  heat  transfer  area  per  chip  is  increased.  With  the  smallest 
channel  size  of  35  pm,  a  minimal  heat  transfer  resistance  of 
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Fig.  4.  Measured  (markers)  and  computed  (lines)  performance  of  pHTSs  with  different 
transfer  resistance  Rhts  and  b)  pumping  power. 


b) 


0.175  cm2I</W  could  be  reached  at  a  fluid  flow  of  0.13  l/min.  The 
measured  Rhts  also  includes  the  resistance  of  the  thermal  paste 
located  between  the  pHTS  and  the  Cu  heater  (approx.  0.08  cm2K/W). 

The  power  necessary  to  pump  the  fluid  through  the  pHTS  is 
depicted  in  Fig.  4b).  The  pumping  power  increases  with  larger  fluid 
flows  and  smaller  microchannel  dimensions.  The  first  is  a  result  of 
an  increased  hydrodynamic  resistance  due  to  the  increased  flow 
velocity.  The  second  can  be  explained  by  the  higher  shear  forces 
(i.e.  pressure  losses)  induced  by  a  larger  surface  to  volume  ratio. 

The  experimental  results  are  in  good  agreement  with  the 
simulation  results.  The  heat  transfer  resistance  (Fig.  4a))  shows  two 
minor  deviations  from  the  computed  values.  First,  the  measured 
values  are  partly  shifted  upwards  or  downwards.  Second,  the  slope 
at  smaller  fluid  flows  does  not  always  correspond  to  the  theory, 
resulting  in  a  slight  overestimation  of  the  thermal  performance. 
The  first  deviation  mainly  results  from  slightly  different  thicknesses 
of  the  thermal  paste  layer,  since  the  exact  thickness  is  very  hard  to 
control.  The  second  deviation  is  likely  related  to  the  observed  small 
air  bubbles  that  can  get  trapped  inside  the  microchannels  at 
smaller  fluid  flows.  This  might  clog  some  channels  and  thus  reduce 
the  heat  convection.  Although  a  correction  factor  was  introduced  to 
account  for  this  phenomenon,  its  impact  is  random  and  can  vary 
between  samples  and  measurements.  The  deviation  of  the  pump¬ 
ing  power  from  the  computed  values  (Fig.  4b))  can  be  explained  by 
the  limited  fabrication  accuracy  of  the  device.  Especially  the  con¬ 
necting  slots  between  the  top  manifold  layer  and  bottom  micro¬ 
channels  are  likely  sources  of  geometric  deviations,  due  to  the 
bonding  and  slot  release  processes  [13]. 

In  order  to  quantify  the  influence  of  the  boundary  conditions  on 
the  thermal  and  hydrodynamic  properties  of  the  pHTS,  the  fluid 
inlet  temperature  and  the  heater  temperature  were  varied  as  well. 
By  increasing  the  fluid  temperature,  an  improvement  of  the  overall 
system  performance  was  observed.  As  a  result  of  an  inlet  temper¬ 
ature  increase  of  2K  (Tin  =  295K),  the  Rhts  measured  at  a  fluid  flow 
of  0.04  l/min  decreases  by  6%  and  the  pumping  power  by  1% 
(averaged  over  3  different  samples).  The  reduction  in  Rhts  origi¬ 
nates  from  the  increase  of  the  thermal  conductivity  of  the  fluid, 
resulting  in  an  increased  thermal  convection  at  the  microchannel 
walls  and  thus  reduced  convective  resistance.  The  higher  inlet 
temperatures  of  the  fluid  also  reduce  the  mean  fluid  density  and 
mean  fluid  kinetic  viscosity.  This  decreases  the  hydrodynamic 


resistance,  leading  to  a  smaller  pressure  drop  and  consequently 
smaller  pumping  powers. 

The  same  trend  is  observed  when  the  heater  temperature  is 
increased.  However,  the  effect  is  less  pronounced,  since  the  fluid  is 
heated  up  indirectly.  In  this  case,  an  average  improvement  of  1%  of 
the  heat  transfer  resistance  and  6%  of  Ppump  was  measured  for  a 
heater  temperature  increase  by  10K  (Thot  =  308K)  and  a  fluid  flow 
of  0.04  l/min  (averaged  over  3  different  samples).  The  measure¬ 
ment  results  demonstrate  an  improved  thermal  and  hydrodynamic 
pHTS  performance  at  higher  operation  temperatures.  However,  in 
the  small  investigated  temperature  ranges,  the  influence  on  the 
pHTS  parameters  is  rather  limited.  Additionally,  due  to  the  mea¬ 
surement  uncertainties,  the  given  changes  should  be  regarded  as 
qualitative  trends. 

3.2.  pTEG/ pHTS  performance 

In  order  to  evaluate  the  optimal  pHTS/pTEG  system  perfor¬ 
mance,  the  influence  of  the  working  conditions  (i.e.  the  applied 
thermal  gradients),  the  microchannel  dimensions  and  the  pTEG 
thickness  (i.e.  thermal  resistance)  on  the  system’s  power  output 
were  investigated. 

3.2.1.  Influence  of  working  conditions 

Fig.  5a)  shows  the  measured  Seebeck  voltage  and  the  output 
power  of  the  pTEG/pHTS  system  with  respect  to  the  applied  tem¬ 
perature  gradient  and  different  fluid  flow  rates.  In  Fig.  5b),  the 
corresponding  net  output  power  is  depicted.  The  slope  of  the 
Seebeck  voltage  (dashed  line)  increases  with  increasing  flow  rates, 
i.e.  decreasing  heat  transfer  resistance  Rhts.  due  to  an  increased 
heat  flux  and  therefore  larger  actual  temperature  gradient  across 
the  generator.  However,  the  differences  in  the  slopes  are  rather 
small  due  to  a  large  pTEG  thermal  resistance  Rteg  compared  to  the 
Rhts  variations,  thus  resulting  in  a  relatively  small  influence  of  the 
enhanced  heat  transfer.  The  specific  influence  of  the  Rteg  to  Rhts 
ratio  will  be  discussed  in  more  detail  in  Section  3.2.3.  Corre¬ 
spondingly,  the  differences  in  the  pTEG  output  power  (solid  lines) 
shows  only  small  deviations  between  different  flow  rates. 

If  the  consumed  pumping  power  is  subtracted  from  the  output 
power  (Fig.  5b)),  the  influence  of  the  flow  rates  increases  signifi¬ 
cantly.  Larger  flows,  consuming  higher  pumping  powers,  need 


750 


V.  Wojtas  et  al.  /  Renewable  Energy  60  (2013)  746—753 


higher  temperature  gradients  to  achieve  a  positive  net  output 
power.  However,  their  slopes  are  steeper  and  therefore,  at  larger 
applied  AT,  higher  flow  rates  outperform  lower  ones.  The  transition 
points  for  this  system  are  marked  in  the  graph  by  vertical  dashed 
lines.  This  means  that  the  system  performance  can  be  adapted  to 
different  operating  conditions  by  varying  the  fluid  flow  rate.  For 
operating  conditions  with  small  temperature  gradients,  small  flow 
rates,  i.e.  large  heat  transfer  resistances  and  small  pumping  powers, 
are  advantageous.  With  increasing  temperature  gradients,  the  flow 
must  be  increased  in  order  to  track  the  maximal  net  output  power 
of  the  system. 

Fig.  6  shows  the  net  output  power  for  pTEGs  with  different  ZT 
values  as  a  function  of  the  applied  temperature  gradient.  At  the 
maximal  applied  temperature  difference  of  95I<,  a  net  output  power 
of  126.3  ±  1.5  mW/cm2  was  measured  with  a  module  ZT  =  0.1  of  the 
pTEG,  resulting  in  an  overall  system  thermoelectric  efficiency  factor 
of  14  pW/K2cm2.  To  show  the  performance  potential  of  devices 
with  improved  thermoelectric  material  properties,  a  computed  net 
power  for  a  module  ZT  =  1  is  added  to  the  graph. 


3.2.2.  Influence  of  microchannel  dimensions 

A  decrease  in  the  microchannel  width  reduces  the  heat  transfer 
resistance  and  increases  the  consumed  pumping  power  of  the  pHTS 
(see  Fig.  4a)  and  b)).  The  influence  of  the  microchannel  dimensions 
on  the  net  output  power  of  the  generator  is  plotted  for  different 
flow  rates  in  Fig.  7.  A  reduction  in  channel  width  directly  leads  to  an 
enhancement  of  the  actual  temperature  gradient  across  the  pTEG 
and  therefore  an  increase  in  the  pTEG's  output  power.  However, 
considering  the  net  pTEG  output  power,  the  simultaneous  increase 
in  required  pumping  power  PPUmp  counteracts  this  positive  effect. 
Therefore,  at  higher  fluid  flows  where  higher  pumping  powers  are 
needed  (blue  line),  Ppump  outweighs  the  reduced  Rhts  for  small 
microchannel  widths,  and  larger  channel  widths  yield  better  net 
power  results.  With  decreasing  flow  rates  the  optimal  channel 
width  shifts  towards  smaller  microchannel  dimensions.  At 
AT  =  30K,  a  fluid  flow  of  0.04  1/min  and  a  channel  width  of  50  pm 
results  in  the  maximal  net  output  power.  Considering  the 
decreasing  relative  influence  of  the  pumping  power  for  increasing 


with  markers)  and  computed  for  a  pTEG  with  ZT  =  1  (solid  line)  as  a  function  of  the 
applied  AT  at  a  fluid  flow  rate  of  0.07  1/min  and  wch  =  35  pm. 


system  as  a  function  of  the  microchannel  width  and  the  applied  fluid  flow  rate 
measured  at  a  AT  of  30K. 
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temperature  differences  AT  applied,  the  impact  of  a  reduced  Rhts 
will  be  more  pronounced.  According  to  simulations,  at  an  applied 
temperature  gradient  of  90K,  the  maximal  net  power  output  is 
reached  with  a  channel  width  of  35  pm,  at  a  volumetric  flow  rate  of 
0.06  1/min.  Therefore,  knowing  the  standard  operating  conditions 
of  the  system,  a  custom  microchannel  design  can  be  applied  to 
achieve  an  optimal  output  performance. 

3.2.3.  Influence  of  fiTEG  thermal  resistance 

As  shown  in  Fig.  5a),  the  influence  of  the  varying  heat  transfer 
resistances,  i.e.  varying  fluid  flow,  on  the  pTEG  output  power  was 
small.  This  can  be  attributed  to  the  fact  that  the  thermal  resistance  of 
the  pTEG  Rteg  =  5.1  cm2I</W)  is  significantly  larger  than  the  variation 
in  the  heat  transfer  resistance  of  the  pHTS  (Rhts  =  0.11-0.48  cm2K/ 
W).  To  investigate  the  influence  of  a  smaller  difference  in  thermal 
resistance,  the  pTEG  thickness  was  reduced  in  two  steps.  The  mea¬ 
surement  results  and  corresponding  simulations  are  summarized  in 
Fig.  8.  The  Seebeck  voltage  decreases  with  decreasing  pTEG  thickness 
and  reduced  flow  rates  (Fig.  8a)).  This  is  a  direct  consequence  of  the 
reduced  temperature  gradient  across  the  generator.  The  Seebeck 
voltage  of  a  generator  can  be  expressed  in  a  simplified  model  (Joule 
heating  and  Peltier  effects  are  neglected)  [28]: 


Vseebeck  =  SATteg  =  S  AT  (1) 

Rteg  +  «con 

where  S  is  the  Seebeck  coefficient  of  the  pTEG,  ATjec  is  the  tem¬ 
perature  gradient  at  the  generator  surface  and  AT  is  the  external 
thermal  gradient.  Rcon  is  defined  as  the  sum  of  all  thermal  contact 
resistances,  including  the  heat  transfer  resistance  of  the  pFITS  Rhts 
and  the  interface  resistances  Rinterface  at  both  sides  of  the  TEG  (see 
Fig.  3 ).  Each  Rinterface  is  estimated  to  be  approximately  0.22  cm2I</W, 
resulting  in  a  total  contact  resistance  Rc0n  between  0.55  cm2I</W 
(0.1 1/min)  and  0.92  cm2K/W  (0.02  1/min). 

At  pTEG  thicknesses  below  50  pm,  the  slope  of  the  Seebeck 
voltage  is  very  steep,  since  the  pTEG  thermal  resistance  is  smaller 
than  the  total  contact  resistance  (Rteg  <  Rcon).  Therefore,  a  small 
increase  in  the  pTEG  thickness  results  in  a  relative  large  increase  in 
ATteg-  With  increasing  thicknesses  (Rteg  >  Rcon),  the  curve  be¬ 
comes  progressively  flat  due  to  the  reduced  influence  of  the  pTEG 
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thickness  on  ATteg-  At  the  same  time,  the  influence  of  varying  Rhts 
becomes  smaller  with  increasing  pTEG  thickness,  which  is  indi¬ 
cated  by  the  smaller  spreading  of  the  Seebeck  voltage  measure¬ 
ments  for  different  fluid  flows.  This  effect  is  more  pronounced  for 
the  pTEG  output  power  (Fig.  8b)),  where  the  Seebeck  voltage  is 
squared  (see  equation  in  Fig.  8b)  inset).  The  pTEG  output  power  also 
increases  with  increasing  pTEG  thickness,  it  reaches  a  maximum 
and  decreases  again.  In  the  first  section,  the  power  increase  is  again 
attributed  to  the  increase  of  ATteg-  At  the  maximum  point,  the 
thermal  resistances  of  pTEG  and  contacts  are  approximately 
matched  [12],  However,  considering  the  Joule  heating  and  Peltier 
cooling  losses,  there  is  a  shift  of  the  optimal  point  to  larger  pTEG 
thicknesses  [28],  On  the  right  side  of  the  maxima,  the  negative 
influence  of  the  increased  electrical  resistance  of  the  generator 
dominates  over  the  increase  in  the  Seebeck  voltage.  With  increased 
fluid  flows  applied,  i.e.  reduced  Rhts.  the  point  of  maximal  output 
power  shifts  towards  higher  values  and  smaller  pTEG  thicknesses.  A 
reduced  heat  transfer  resistance  results  in  higher  ATTeg  and  smaller 
pTEG  thicknesses  to  achieve  thermally  matched  resistances. 

For  each  temperature  operating  range  of  the  system,  an  optimal 
pTEG  thickness  can  be  defined  to  maximize  the  power  output.  For 
applications  with  lower  AT  available,  where  smaller  fluid  flows  can 
be  applied  to  reach  the  maximum  net  output  power  (see  Fig.  8b)), 
the  optimal  pTEG  thickness  is  around  90  pm.  With  increasing 
external  temperature  gradient  and  therefore  increasing  required 
fluid  flows  for  optimal  power  exploitation,  the  pTEG  thickness  can 
be  reduced  even  further.  The  presented  optimal  pTEG  thicknesses 
are  specific  values  for  the  pTEG  used  in  this  study.  Changes  in 
thermal  and  electrical  properties  of  the  generator  will  result  in  a 
shift  of  the  maximal  power  points  to  different  pTEG  thicknesses. 
Small  thicknesses  enabled  by  the  very  low  heat  transfer  resistance 
pHTS  are  of  particular  interest  for  energy  harvesting  from  waste 
heat  recovery.  For  example  in  industrial  heat  exchange  systems, 
where  high  heat  fluxes  are  required  to  sustain  the  necessary  heat 
dissipation,  a  low  thermal  resistance  of  the  generator  is  favorable. 
Additionally,  a  reduced  pTEG  thickness  brings  the  advantages  of 
smaller  material  and  fabrication  costs,  as  well  reductions  in  system 
size  and  weight. 

Another  important  aspect  related  to  low  heat  transfer  re¬ 
sistances  is  the  relevance  of  the  thermal  properties  of  the  ther¬ 
moelectric  material.  The  thermal  resistance  of  a  generator  depends 


d  the  thermal  contact  resistance  Ro„  (0-5  cm2K/Wm,  in 
:he  value  of  X  =  0.2  W/mK.  A  pTEG  thickness  of  215  pm, 
ty  of  the  passive  material  of  0.3  W/mK,  and  parameters 


not  only  on  the  thickness,  but  also  on  the  thermal  conductivity  of 
the  active  (and  passive)  material.  When  working  with  very  small 
thermal  contact  resistances,  the  importance  of  low  thermal  con¬ 
ductivities  is  reduced.  In  Fig.  9  the  relative  pTEG  output  power  is 
plotted  versus  the  thermal  conductivity  A  of  the  thermoelectric 
material,  for  different  thermal  contact  resistances  Rcon-  The  relative 
power  is  normalized  to  a  pTEG  output  power  value  potentially 
achieved  with  a  very  low  thermal  conductivity  of  A  =  0.2  W/mK 
(e.g.  of  a  superlattice  structure  [8]).  The  relative  influence  of  an 
increase  in  A  decreases  with  reduced  Rcon-  In  the  ideal  case  of 
i?con  -»■  0,  the  output  power  becomes  a  function  of  only  the  Seebeck 
coefficient  S,  the  electrical  resistance  Rei  and  the  externally  applied 
AT  (see  Formula  1  and  inset  in  Fig.  8b)).  Therefore,  when  working 
with  very  low  heat  transfer  resistances,  the  material  optimization 
efforts  should  shift  from  the  optimization  of  the  figure  of  merit 
ZT  —  u2(jTI'a  to  the  enhancement  of  the  power  factor  P  —  a2tr.  This 
would  in  particular  refer  to  nanostructured  composites  [29,30]  and 
superlattice  structures  [8],  where  an  improvement  of  ZT  is  mainly 
achieved  by  the  reduction  of  A.  This  conclusion  is  in  agreement  with 
Narducci  [31  ],  who  also  emphasized  the  much  higher  relevance  of 
the  power  factors  over  the  ZT  value,  especially  for  high  heat  flux 
applications. 


pumping  power  outweighs  the  reduced  Rhts  and  thus  smaller  flows 
and  wider  channels  lead  to  higher  net  powers.  With  a  larger  tem¬ 
perature  gradient  applied,  the  influence  of  the  consumed  pumping 
power  decreases  and  the  optimal  flow  rate  increases,  whereas  the 
optimal  microchannel  width  decreases. 

A  system  net  output  power  of  126.3  ±  1.5  mW/cm2  was 
measured  with  a  module  ZT  of  0.1  at  the  applied  temperature  dif¬ 
ference  of  95K  and  a  fluid  flow  rate  of  0.07  1/min.  It  was  shown  that 
with  the  achieved  thermal  coupling,  the  pTEG  thickness  could  be 
further  reduced  to  60—90  pm  in  order  maximize  the  output  power 
under  thermally  matched  resistances.  This  brings  several  advan¬ 
tages  with  respect  to  fabrication  and  cost  and  makes  the  system 
suitable  for  high  heat  flux  waste  heat  recovery  applications. 
Moreover,  the  reduced  relevance  of  the  thermal  conductivity  in  the 
thermoelectric  material  with  decreasing  contact  resistances  was 
highlighted.  The  discussed  advantages  and  achieved  performance 
demonstrated  the  potential  for  future  up-scaling  and  further  sys¬ 
tem  extension  into  compact  multi-layer  heat  exchangers  for  in¬ 
dustrial  applications. 
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4.  Conclusions 


The  proposed  system  integration  of  pTEGs  with  a  low  power 
consuming  pHTS  was  characterized  with  respect  to  the  generated 
net  output  power.  The  small  heat  transfer  resistances  enabled  much 
higher  effective  temperature  gradients  across  the  pTEGs  and  led  to 
a  significant  enhancement  of  the  produced  output  power.  The 
presented  pHTS  achieved  heat  transfer  resistances  between  0.17 
and  0.88  cm2K/W  at  corresponding  low  pumping  powers  of  6.8  - 
0.06  mW/cm2. 

Due  to  the  coupled  influence  of  pHTS/pTEG  system  parameters 
and  boundary  conditions,  its  net  output  performance  can  be  opti¬ 
mized  for  different  operating  conditions  by  the  applied  fluid  flow 
rate  and  by  the  careful  design  of  the  microchannels.  A  trade-off 
must  be  made  between  a  low  heat  transfer  resistance,  i.e.  a  high 
effective  temperature  gradient,  and  the  consumed  pumping  power. 
At  smaller  operating  temperature  differences,  the  impact  of  the 


